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Abstract. We present a new superfast algorithm for solving Toeplitz systems. This algorithm is 
based on a relation between the solution of such problems and syzygies of polynomials or moving 
CO , lines. We show an explicit connection between the generators of a Toeplitz matrix and the generators 

of the corresponding module of syzygies. We show that this module is generated by two elements and 
' the solution of a Toeplitz system T u = g can be reinterpreted as the remainder of a vector depending 

s ! ' on g, by these two generators. We obtain these generators and this remainder with computational 

£H i complexity C(nlog 2 n) for a Toeplitz matrix of size n X n. 

a 
s 

CN 1. Introduction 

Structured matrices appear in various domains, such as scientific computing, signal processing, 
. . . They usually express, in a linearize way, a problem which depends on fewer parameters than the 
number of entries of the corresponding matrix. An important area of research is devoted to the 
development of methods for the treatment of such matrices, which depend on the parameters defining 
them. 

Among well-known structured matrices, Toeplitz and Hankel structures have been intensively stud- 
ied [H H2] . Nearly optimal algorithms are known for their multiplication by a vector and the solution of 
linear systems, for such structure. Namely, if A is a Toeplitz matrix of size n, multiplying it by a vector 
or solving a linear system with A requires 0(n) arithmetic operations (where 0(n) = 0(n\og c (n)) for 
some c > 0) [21 [SO] - Such algorithms are called superfast, in opposition with fast algorithms requiring 
if} \ 0(n 2 ) arithmetic operations. 

The fundamental ingredients in these algorithms are the so-called generators [12] . encoding the 
minimal information stored in these matrices, and on which the matrix transformations are translated. 
The correlation with other types of structured matrices has also been well developed in the literature 
|171 116) . allowing to treat efficiently other structures such as Vandermonde or Cauchy-like structures. 

Such problems are strongly connected to polynomial problems [5j[T]. For instance, the product of 
a Toeplitz matrix by a vector can be deduced from the product of two univariate polynomials, and 
thus can be computed efficiently by evaluation-interpolation techniques, based on FFT. The inverse 
of a Hankel or Toeplitz matrix is connected to the Bczoutian of the polynomials associated to their 
generators. Such a construction is related to Gohberg-Semencul formula [7] (or Trench algorithm 
|18j). which describes the inverse of a Toeplitz matrix in terms of the solution of two specific Toeplitz 
systems (see also Gohberg-Prupnick formula [6]). 

Most of these methods involve univariate polynomials. So far, few investigations have been pursued 
for the treatment of multilevel structured matrices [191 113) . related to multivariate problems. Such 
linear systems appear for instance in resultant or in residue constructions, in normal form computa- 
tions, or more generally in multivariate polynomial algebra. We refer to 15J for a general description 
of multi-structured matrices and their correlations with multivariate polynomials. Surprisingly, these 
multivariate structure also appear in numerical scheme and preconditionners [13] . A main challenge 
here is to devise superfast algorithms of complexity 0(n) for the solution of multi-structured systems 
of size n. 

In this paper, we re-investigate the solution of Toeplitz systems T u = g from a new point of view 
which can be generalized to two-level Toeplitz systems. We correlate the solution of such problems with 
syzygies of polynomials. We show an explicit connection between the generators of a Toeplitz matrix 
and the generators of the corresponding module of syzygies. We show that this module is generated 
by two elements of degree n and the solution of T u = g can be reinterpreted as the remainder of an 
explicit polynomial vector depending on g, by these two generators. We give two algorithms, with 
computational complexity 0(rolog 2 n), to compute the generators of the module of syzygies. We give 
finally an algorithm, with computational complexity 0(n\og n), for the division of the generators 
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by the polynomial vector depending on g. Our new syzygy approach can be connected with Pade 
approximation method developed in [3] to compute efficiently particular solutions of Toeplitz linear 
system. But we replace the computation of generators of structured matrices by the computation 
of generators of a syzygy module and the solution of the linear system from particular solutions by 
Euclidean reduction by the generators of the syzygy module. 

Let R = K[x]. For n £ N, we denote by K[at] n the vector space of polynomials of degree < n. 
Let L = K[x,a; _1 ] be the set of Laurent polynomials in the variable x. For any polynomial p — 
2i=-m Pi x% £ L, we denote by p + the sum of terms with non-negative exponents: p + = Y^7=o Pi x% an< ^ 
by p~ , the sum of terms with strictly negative exponents: p~ = Yl7—-m Pi x% • ^ e nave P = P + + P~ ■ 

For n G N, we denote by il„ = {ui; ui n — 1} the set of roots of unity of order n. 

For a vector u = (uq, . . . , Ufc-i) T S K k , we denote by u(x) the polynomial of degree k — 1 given by 
u(x) = X^iLo 1 u i x% ■ Conversely, if v(x) = ^2^ = q V{X l is a polynomial of degree k — 1, we denote by v 
the vector of length k of coefficients of v(x). 

If no confusion arises, we may also use v to denote the polynomial v{x). 

2. Sygygies and Toeplitz matrices 
Let T G K" x ™ be an n X n Toeplitz matrix. Then T is of the following form: 

I to i-l . . . ^-n+l\ 

t\ to ' ■ : 
: '■■ '•■ t-i 

\tn-l ■•• tl to ) 

Let g = (go, . . . , <? n _i) G IK™ be a vector of length n. We are interested in the following problem: 
Problem 2.1. Find u = (u , . . . , u n -i) G swc/i i/iaf 
(2) T U = g . 

Definition 2.2. Let E = {1, . . . ,cc™ -1 }, anrf LT^ be the projection of L on the vector space generated 
by E, along (x n , x n+1 , . . .) . 

Definition 2.3. From the matrix T and the vectors g and u we define the following polynomials: 



(1) 



T(x)= Yj 

i=— n+1 

2n-l |- 

T(x) = iiX* Witt = < 

n— 1 n—1 



if i < n 
*i-2n ifi>n 



i=0 i=0 

Notice that T(x) is a Laurent polynomial and that T(a;) is a polynomial of degree 2n — 1. By 
construction, we have the following properties: 

Proposition 2.4. f = T+ +x 2n T- andT(w) =f(w) ifwett 2n . 

Proof. We can deduce directly, from the definition of T(x) and T(x), that T = T + +x 2n T~ . Moreover, 
since w 2n = 1 and f(x) = T+(x) + x 2n T~(x), then f(w) = T+(w) + T~(w) = T(w). □ 

According to Proposition 2.1.2 of [TS], we have the following relation between the problem 12. II and 
polynomials: 

Proposition 2.5. We have 

T 'u = g He(T(x)u(x)) = g(x). 

As He(T(x)u(x)) is the polynomial T(x)u{x) from which we remove terms of negative degree and 
of degree > n, then we can write T(x) u(x) as following: 

Proposition 2.6. 

(3) T(x) u{x) = U E (T(x)u(x)) + x- n A(x) + x n B(x), 

where A(x) G K[a;]„-i and B(x) G K[a;] n -2- 
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Proof. By expanding T(x)u(x) we can write 

T(x)u(x) = U E (T(x)u(x)) + (a_ n+1 x- n+1 + ■■■ + a^aT 1 ) + (a n x n + ■■■ + a 2n _ 2 x 2n ~ 2 ) 
= U E (T(x)u(x)) + x~ n (a- n+1 x + ■■■ + a^x^ 1 ) + x n (a n + ■■■ + a 2n _ 2 x n ~ 2 ) 
= U E (T(x)u(x)) + x- n A(x) + x n B(x) 

□ 

Therefore, according to Proposition 12.51 and Proposition 12.61 if u is solution of Tu — g then there 
exist two polynomials A(x) and B(x) in K[a;] n _i such that 

(4) T{x)u{x) - x- n A{x) - x n B{x) = g{x). 

By evaluation at the roots w G ii 2n , and since u~ n — uj n and T(lo) — T(uj) for uj € il 2n , we have 

T(uj)u(uj) + ui n v(ui) — g(ui),Vu) G il 2n (oj), 

where v(x) = —A(x) — B(x) of degree < n — 1. Therefore the polynomial T(x)u(x) + x n v(x) — g(x) 
is multiple of x 2n — 1. We deduce that there exists w(x) G K[x] such that 

(5) f{x)u(x) + x n v{x) + {x 2n - l)w(x) = g(x). 

Notice that w(x) is of degree < n — 1 because (x 2n — 1) w(x) is of degree < 3n — 1. 

2.1. Syzygies. The solutions of Equation (JSJ is a particular case of the following problem, related to 
interesting questions in Effective Algebraic Geometry. 

Problem 2.7. Given three polynomials a,b,c G R respectively of degree < l,< to, < n, find three 
polynomials p,q,r G R of degree < v — l,< v — m,< v — n, such that 

(6) a(x)p(x) + b{x) q(x) + c(x) r(x) — 0. 

The polynomial vector (p, q, r) G K[x] 3 is called a syzygy of (a, b, c). We denote by C(a, b, c) the set 
of syzygies (p,q,r) G K[x] 3 of (a,b,c), i.e. the solutions of ([6]). It is a K[x]-module of K[x] 3 and it 
is called the module of syzygies of (a, b, c). The solutions of Problem 12.71 are £(a, b, c) n K^]^;-! x 
K[x] v - m -i x K[a;]„_ n _i. 

Given a new polynomial d(x) G K[x], we denote by C(a, b, c; d) the set of (p, q, r) G K[x] 3 such that 

(7) a(x)p(x) + b(x) q{x) + c(x) r(x) = d(x). 

Theorem 2.8. For any non-zero vector of polynomials (a, 6, c) G K[a;] 3 , the ¥L[x]-module C(a,b,c) is 
free of rank 2 . 

Proof. By the Hilbert's theorem, the ideal / generated by (a,b,c) has a free resolution of length at 
most 1 (see [H chap. 6] ) , that is of the form: 

->■ K[x] p K[x} 3 -» K[x\ -)• K[x]/I -)• 0. 

As / ^ 0, for dimensional reasons, we must have p — 3 + 1 = 0, then p = 2. □ 

Definition 2.9. For a polynomial vector p — (pi, . . . ,pk) G K[ir] fe , we define 

deg(pi ,...,Pk) = max(deg(pi),...,deg(p fe )). 

Definition 2.10. Assume that deg(p, q, r) < deg(p', q' , r'). A fi-base of C(a,b,c) is a basis {(p,q,r), (p',q', 
of C(a, b, c), with deg(p, q, r) = /i. 

We have the following relation between the degrees of the two elements of a basis of C(a, b, c): 

Proposition 2.11. Let {(pi, qi, Ti), (p 2 , q 2 , r 2 )} be a basis of C(a, b, c), [i\ = deg(pi, qi,r±) and \x 2 = 
deg(p 2 ,q 2 ,r 2 ). We have deg(a,6, c) = + [i 2 . 

Proof. We have 

->• Klx^-d-^ © Klxlu-a-^ ->• K[x\l_ d -)• K[x] v K[x] u /(a, b, c) v -> 0, 

for v ^> 0. As the alternate sum of the dimension of the K-vector spaces is zero and K[x)^/(a,b,c) l , 
is for v ^S> 0, we have 

= 3 (d — v — I) + v — Hi — d + 1 + v — /j,2 — d+l + v+ l = d — fii — \i 2 . 

□ 
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2.2. The module C(f(x),x n , x 2n — 1). Returning to the initial problem, we saw that if u is solution 
of Tu — g then there exist two polynomials v(x) and w(x) in K[x] n _i such that (u(x), v(x), w(x)) 
G C(f(x),x n ,x 2n - l;g(x)). 

By the proposition ^. 11[ if (p, q, r) and(p', q', r') form a basis of C(T(x), x n ,x 2n — l) of degree ui and 
\i 2 respectively then we have ui + u 2 = 2n. We are going to show now that in fact C(T(x), x n ,x 2n — 1) 
has a n-basis, that is a basis of two elements of degree «i = /i 2 = n: 

Proposition 2.12. The K[x]-module C(T (x) , x n , x 2n — 1) has a n-basis. 
Qfpof. Consider the linear mqg[ x ]3_ i K[x] 3 „_i 

(p(ar),g(ar),r(x)) i-> f (x)p(ar) + a; n g(x) + (x 2n - l)r(x), 
which 3n x 3n matrix is of the form 

(9) S := 

where T ,Ti,T 2 are the coefficient matrices of (T(x), xT(x), . . . , x n T(x)), respectively for the list of 
monomials (1, . . . , x™ _1 ), (x n , . . . , x 2 "" 1 ), (a; 2 ™, . . . , a; 3 ™" 1 ). Notice in particular that T = T + T 2 . 

Reducing the first block (To|0| — I n ) by the last block (T2 1 1 1„ ) , we replace it by the block (To + 
T2|0|0), without changing the rank of S. As T = To + T2 is invertible, this shows that the matrix S 
is of rank 3n. Therefore ker(5 l ) = and there is no syzygies in degree n — 1. 

As the sum 2n = [ii+[i 2 , where fi 2 are the degrees of a pair of generators of £(T(x), x n ,x 2n — 1), 
and as [i\ > n and [i 2 > n, we have fxi = fi 2 = n. Moreover, £(T(x), x n , x 2n — 1) is free of rank 2. 
Thus there exist two linearly independent syzygies {u\, v\,w\), (u 2l v 2 ,w 2 ) of degree n, which generate 
C(f (x),x n ,x 2n - 1). □ 

A similar result can also be found in [20j . but the proof much longer than this one, is based on 
interpolation techniques and explicit computations. 

Let us now describe how to construct explicitly two generators (iti, vi,wi), {u 2l v 2 ,w 2 ) of C(T(x), x n , x 2n - 
1) of degree n. As T(x) is of degree < 2n — 1 and the map §8§ is surjective, there exists (u,v,w) G 
K^ai] 3 ^! such that 

(10) f (x)u(x) + x n v{x) + (x 2n - l)w = f(x)x n . 

We deduce that (ui,Vi,wi) = (x n — u, —v, -w) G C(T(x), x n ,x 2n - 1). 
Since there exists (u', v' , w') G Kfa;]^! such that 

(11) f(x)u'(x) + x n v'(x) + (x 2n -l)w' = l = x n x n - (x 2n - 1), 

we deduce that (w 2 , v 2 , w 2 ) — (— u', x n - v', —w' — 1) G C(T(x), x n , x 2n - 1). 

Now, T and are linearly independent since by construction, The coefficient vectors of x" in 
(ui, Vt, Wi) and (u 2 , v 2 , w 2 ) are respectively (1, 0, 0) and (0, 1, 0), which shows that vectors (iii,Ui,u>i), 
(u 2 ,v 2 ,w 2 ) G C(T(x), x n , x 2n — 1) fl K[i]„ are linearly independent. Therefore, they form a basis of 
C(f (x),x n ,x 2n - 1). 

Now we can prove our aim theorem: 

Theorem 2.13. The vector u is solution of ^ if and only if there exist v(x) and w(x) in K[x] n _i 
such that 

(u(x),v(x),w(x)) G C(f {x),x n ,x 2n - l;g(x)) 

Proof. If u is solution of (J^J), we see that there exist v(x) G JK[x]„_i and w{x) G K[x]„_i such that 

f{x)u{x) + x n v{x) + {x 2n - l)w(x) = g(x). 

Conversely, a solution (u(x), v(x), w(x)) G C(T(x), x n , x 2n — 1; g(x))tlK[x]f l _ 1 implies that (u, v, w) G 
K 3n is a solution of the linear system: 

S 

where S is has the block structure ©, so that T2 u + w = and To u — w = (To + T 2 )u = g. As we 
have To + T 2 = T, the vector u is a solution of ([2]), which ends the proof of the theorem. □ 
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Computing the inverse of a Toeplitz matrix T is equivalent to computing the first and the last 
column of T -1 , based on Gohberg-Semencul decomposition (see [51 HH UHl E] for more details about 
Gohberg-Semencul decomposition) . 

We are going to show that the solutions of Equations (fit))) and (fTTj) which gives us the n-basis 
{(ui,Vi,Wi), (u2,V2,W2)} is related to the solution of two specific Toeplitz linear systems. 

Proposition 2.14. Let (u(x),v(x),w(x)) and (u'(x),v'(x),w'(x)) be mK n _i[x] such that 

f f(x)u(x) + x n v(x) + (x 2n - 1) w(x) = f(x)x n , 
\ f(x)u'(x) + x n v'(x) + (x 2n - 1) w'(x) = 1. 

Then Tu' — e\ and Tu — ZTe n , with Z is the lower shift matrix of size n. 

Proof. As u'(x), v'(x), w'(x) and 1 are of degree < n— 1, then, by Theorem l2.131 T(x)u'(x)+x n v'(x) + 
(x 2n — 1) w'(x) = 1 is equivalent to Tu' = e\ (ex(x) — 1) and u' is the first column of T _1 . 
We have f{x) = T+(x) + x 2n T_(x), then 

f(x)u(x) + x n v(x) + (x 2n - l)w(x) = x n T+(x) + x n ((x 2n - l)T_(a) + T_(x)). 

Therefore, 

f (x)u(x) + x n (v(x) - T+(x)) + (x 2n - l)(w{x) - x n T_(x)) = x n T_(x). 

As x n T_(x) is of degree < n — 1 and is the polynomial associated with the vector ZTe n , by Theorem 
u is such that Tu = ZTe n . □ 

Notice that u is not the last column of T _1 , but we can use u and v! to compute it (see [5]). 
Therefore, defining a n-basis of C(T (x) , x n , x 2n — 1) from the solution of Equations (TTu| and pT|) is 
equivalent to computing the Gohberg-Semencul decomposition of T _1 . 

In the following section, we reduce translation of the solution of Tu — g to an Euclidean divi- 
sion, based on our decomposition, instead of multiplying g by triangular Toeplitz matrices, based 
on Gohberg-Semencul decomposition. The advantage of our decomposition is that we can general- 
ized it to two-level problems, which allows us to describe a "Gohberg-Semencul" decomposition of 
Toeplitz-block- Toeplitz matrices. 



3. Euclidean division 

In this section, we show how to obtain the solution vector (u(x),v(x),w(x)) G C(T(x),x n ,x 2n — 
1; g{xj) n K[x]^ from a n-basis of £(T(x), x n , x 2n — 1) and a particular solution in C(T(x), x n , x 2n — 

!;<?(*))• 

From Theorem 12. 131 we deduce the two following corollaries: 

Corollary 3.1. For all g{x) £ K n _i[x], the set £,(T(x),x n 7 x 2n — l;g(x)) ("1 K^_ 1 [ac] has exactly one 
element. 

Proof. As T is invertible, there exists a unique u such that Tu = g. From the theorem I2.13[ there 
exists v(x), w(x) of degree < n— 1, such that (u(x), v(x),w(x)) £ jC(T(x), x n ,x 2n — 1; g(x))nK'} l _ 1 [x}. 

The uniqueness is also obvious: if (u'(x), v'(x), w'{x)) € £(T(x), x n , x 2n — l;g(xj) flKj.Ji], then 
(u(x),v(x),w(x)) - (u'(x),v'(x),w'(x)) G C(f(x),x n ,x 2n - 1) n^fi] which equal {(0, 0, 0)} (see 
the demonstration of the proposition ^. 121) . Then (u(x), v(x), w(x)) = (u'(x),v'(x),w'(x)). □ 

Corollary 3.2. Let {(ui,vi,w{), (u2, 1^2,102)} be a n-basis of £(T(x),x n ,x 2n — 1). Let{p,q,r) be in 
C{T{x), x n , x 2n — l;g(x)). There exists a unique (u,v,w) € L(T(x) 1 x n ,x 2n — l;g(x)) PiK.^_ 1 [x] and 
a unique pair of polynomials p\ and P2 such that 

>r\ M M ( 
q = pi \ v i \ + P2 \ v 2 + \v 

' \ w lJ \ W 2/ \ W J 

This decomposition is called the division 0/ (p, q,r) by (ui,Vi,Wi) and (112, V2, 102). 

Proof. From the previous corollary, there exist a unique element in C(T(x), x n , x 2n — 1; g(x))nK.f l _ 1 [x], 
let (u, v, w) be this element. As {(mi, v±, Wi), {U2, V2, 102)} is a n-basis of C{T{x), x n , x 2n — 1), and as 



6 HOUSSAM KHALIL, BERNARD MOURRAIN, AND MICHELLE SCHATZMAN 

(p, q, r) — [u, v, w) s C(T(x), x n , x 2n — 1), then there exist a unique pair of polynomials unique p\ and 
p 2 such that 

' m \ i u 2 

Pi I «1 +P2 "2 



VU>2, 



□ 



As a consequence of the two corollaries, we have the following important property: 

Theorem 3.3. Let {(u±, v\, Wi), (tt2, V2, W2)} be an-basis of £(T(x),x n ,x 2n — l), and let g £ 1" . The 

I \ /m, u 2 \ 

remainder of the division of \ x n g I by \ Vi v 2 is the unique element (u,v,w) £ C(T (x) , x n , x 2n — 

\ 9 J w 2 ) 
\;g(x)) flK^Ji], and therefore u is the solution of Tu = g. 



Proof. The vector I or™ <? S £(T(x), x n , x 2n — l;g) is a particular solution. We reduce it by 



-9 , 



U\ u 2 

vi v 2 I and obtain 
,u>i w 2 . 



\ / iti u 2 
g - «i w 2 
3 / \wi w 2 . 



where (u,v,w) € Kfa;]^! n C{T(x), x n , x 2n — l;g) is the remainder of division. Thus (u, v, w) is the 
unique vector G K[x]3_ 1 n £(T(a;), x n , x 2n - 1; 5 ). □ 

A way to perform the division is to choose a n-basis {(ui, Vi,u>i), (u 2 ,v 2 , w 2 )} of C(T(x), x n , x 2n — l) 
so that the 2x2 coefficient matrix of x n in 

u\{x) u 2 {x) 
vi(x) v 2 (x) 

is invertible. In this case we can reduce the polynomial (0, x n g(x)) to reach to a degree < n — 1 and 
we can write in a unique way 

By the uniqueness of the remainder in the Euclidean division, we obtain the following proposition: 

Proposition 3.4. The first coordinate of the remainder in the division of ( ? | by ( U2 \ is the 

\x n gj y \vi v 2 J 

polynomial u(x) such that its associated vector u is the solution of T u = g. 
So we set the following problem: 

Problem 3.5. Given a matrix and a vector of polynomials ^j^j f'o) °^ ^ e ^ ree n suc ^ ^at 

e ? e fi \ is invertible and \ ^^\ \ of deqree m > n, find the remainder of the division of \/ X T\ by 
e(x) e'(x) 

m f\x) 

We describe here a generalized Euclidean division algorithm to solve problem 
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Let E(x) = ^(xjl of degree to, B(x) = (j^Jj of de § ree n < m. E(x) = B(x)Q(x)+R(x) 

with deg(R(x)) < n, and deg(Q(x)) < m — n. Let z = -. We have 

E{x) = B(x)Q(x) + R(x) 

* =B( i )Q( I )+i?( i } 

<^ = z n B(-)z ro - n g(-) + z m - n+l z n - 1 R(-) 

z z z z 

(12) <=> E{z) =B(z)Q(z) + z m ~ n+1 R(z) 

with E (z) , B (z) , Q (z) , R(z) are the polynomials obtained by reversing the order of coefficients of 
polynomials E(z), B(z),Q(z), R(z). 

m => B{zf 1 E(z)^Q(z) + z m+n - 1 B(zf 1 R{z) 



m— n+1 



=> Q{z) = B(z) E(z) mod z 

The formal power series B(z) exists because the constant coefficient of B(z) is invertible. Thus 
Q(z) is obtained by computing the first m — n + 1 coefficients of B(z) E(z), which is obtained by 
computing W(x) = B{z) , then by multiplying W(x) by E(z). 

To find W(x) = B(z)^ 1 we use Newton's iteration. Let f(W) = B - W' 1 . We have 

f'{Wi).(Wi +1 - Wi) = -WrHWi+i ~ W^Wf 1 = f(W t ) = B — Wf 1 . 

Thus we set 

W i+ i=2Wi-WiBW l , 
and Wo = Bq , which exists. Moreover, we have 

W - Wi+i = W-2Wi + WiBWt 
= Wih-BWif 
= (W - Wi)B(W - Wi) 
Thus Wi(x) = W(x) mod x 11 for I = 0, . . . , |~log(m - n + 1)1 . 

Proposition 3.6. We need O(n\og(n) log(ro — n) + rn log to) operations to solve vroblem \3.5\ 

Proof. We must do [log(TO — n + 1)] Newton's iteration to obtain the first m — n + 1 coefficients of 
B(z) = W(x). And each iteration requires O(nlogn) operations (multiplication and summation of 
polynomials of degree n). Finally, multiplication B(z) E(z) requires (Dijn log to) operations. □ 

Notice that, for our problem m = n and this algorithm requires 0(nlog 2 n) arithmetic operations. 
In the following section, we show how to compute a n-basis in 0(n log 2 n) arithmetic operations. 

4. Construction of the generators 

The canonical basis of K[x] 3 is denoted by a±, cr%, 03. Let p±, P2 be the generators of £(T(x), x n , x 2n — 
1) of degree n given by 

, , pi = x n ai - (u, v, w) = (ui,vx,wi) 

Pi = x n (7 2 - [u', v', w') = (u 2 , v 2 , W2), 

where (u,w,w), {v! , v' , w') are the vectors given in (flu]) and (fTTj) . 

We describe two methods for computing (ui,ui,tui) and (^,2,^2,^2)- The first one uses the Eu- 
clidean gcd algorithm, the second one is based on the method in [213] . 

We recall firstly the algebraic and computational properties of the well known extended euclidean 
algorithm (see 2 lj ) : Given p(x) , p 1 (x) two polynomials in degree to and to' respectively, let 

ro =P, n= p', 

s = 1, ai = 0, 

< = 0, h = l. 

and define 
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Si+1 



<*i — 1 yi^i; 
U-l — Qiti, 



where qi results when the division algorithm is applied to and rj, i.e. r.j_i = g^r.; + r^ + i . 
Proposition 4.1. Lei i G N sitc/i i/iai r; = 0. Then r;_i = gcd(p(x),p'(x)). 

And more generally we have: 
Proposition 4.2. For all i — 1, ... ,1 we have 

SiP + Up' = n and (si,U) = 1, 

and 

( degr 4+ i < degr-i, i = l,...,l-l 

degSi+i > degSi and degi i+ i > degU, 
degs i+ i = deg(qi.Si) = degv - degn, 
degt i+1 = deg(%.tj) = degu - degrj. 

We can now present our algorithm. It can be found in the proof of the following theorem: 

Theorem 4.3. By applying the Euclidean gcd algorithm to p(x) = x n T andp'(x) = x 2 ™" 1 stopping 
in degree n — 1 and n — 2, we obtain p\ and pi respectively. 

Proof. We see that Tu = g if and only if there exist a(x) and o(x) in K[x] n _i such that 

f(x)u(x) + x 2 "" 1 ^) = x n - l g{x) + a(x), 

where T(x) = x n ~ 1 T(x) is a polynomial of degree < 2n — 2. In (|10p and (fTTj) we saw that for g(x) = 1 
(.9 = e i) an d 5(2^) = x n T(x) (g — (0, t- n +i, . . . , t-i) T ) we obtain a base of C(T(x), x n , x 2 " — 1). 
Notice that Tu\ = e\ if and only if there exist cti(x) G K[x] n _2, &i(x) G K[x]„_i such that 



(14) 



T{x) Ul (x) +x 2n - 1 b 1 (x) 



oi(x), 



and T«2 = (0,t_ n +i, . . . , t_i) T if and only if there exist 02 (x) G K[x]„_2, &2(x) G K[x]„_i such that 
(15) f{x)(u 2 (x) + x n ) + x 2n - x b 2 {x) = a 2 {x). 

As degai(x) < n — 2 and dega2(x) < n — 2, by applying the extended Euclidean algorithm in 
p(x) — x n ~ x T and p'(x) = x 2n_1 until we have degn(x) = n — 1 and degr;+i(x) = n — 2 we obtain 



and 



Ztl(x) = — S;(x), 6i(x) = — <;(x), X™ 1 + di (x) = — n (x) , 
Ci Ci Ci 

' l + u 2 (x) = — s;+i(x), & 2 (x) = — t;+i(x), a 2 (x) = — n+i(x), 

C2 C 2 C 2 



with ci and C2 are the highest coefficients of r;(x) and s; + i(x) respectively. In fact, Equation (jlip is 
equivalent to 



n-l 



n-1 | 


/ i-n+l 






to 








n / 












to 






n-l | 


V 


tn-1 


1 


1 



/ \ 

Ui 

h 

v y 



Ol 



W 



since T is invertible then the (2n — 1) X (2n — 1) block at the bottom is invertible and then u\ and b\ 
are unique. Therefore a\ is also unique. 
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As dcgr; = n— 1 then, by Proposition 14.21 degs/+i = (2n — 1) — (n — 1) = n and degi;+i = 
(2n — 2) — (n — 1) = n — 1. By the same proposition, we also have deg s; < rt — 1 and degi/ < n — 2. 

Therefore, degui = degs; and deg&i = dcgi/. Then as iti(x) and ^-s; are unitaries, ^-s;(a;) = 
Ui(x) which implies that -^ti(x) = bi(x). For the same reasons, we have x n + U2{x) = -^si + i(x) and 
b2{x) = ±ti + i{x). 

Finally, Tu = e\ if and only if there exist v{x), w{x) such that 

(16) f(x)u(x) + x"v(x) + (x 2n - l)w(x) = 1. 
As f(x) = T+ + x 2n T- = T + (x 2n - 1)T~, we deduce that 

(17) T(x)u(x) + x n v(x) + (x 2n - l)(w(x) + T~(x)u(x)) = 1. 

Moreover, we also have T(x)u(x) — x~ n+1 ai(x) +x n bi(x) = 1 and x~ n+ ai(x) = x n (x ai) — x~ n (x 2n — 
l)xai. Thus 

(18) T(x)u(x) + x n (b(x) -xa(x)) + (x 2n - l)x- n+1 a{x) = 1. 

Comparing (IT71) and ([T5]) . and as 1 = x n x n — (x 2n — 1) we deduce that w(x) — x~ n+1 a(x) — T-(x)u(x) + 
1, which is the part of positive degree of — T-(x)u{x) + 1. This conclude the proof of the proposition. 

□ 

Remark 4.4. The usual Euclidean gcd algorithms are of computational complexity 0{n 2 ), but su- 
perfast euclidean gcd algorithms use no more then 0{nlog 2 n) operations, exist. See for example |21j 
chapter 11. 

The second method for computing (tti,«i,tui) and (u2,V2,W2) * s °^ polynomials an d interpolation 
points. We are interested in computing the coefficients of the canonical basis element <n, oi in this 
basis. The coefficients of 03 can be obtained by reduction of (T(x) x n ) B(x) by x 2n — 1 where 

T3( T \ — f U ( X ) U '( X )\ 

{, ~\v{x) v'ix)) 1 

where (u,v), (u',v') are the two first coordinates of the solution of Equations ([TO]) and (fTTjl . A 
superfast algorithm for computing B{x) is given in [2D]. Let us describe how to compute it. 

By evaluation of (fTB]) at the roots uij G i^n, we deduce that (u(x),v(x)) and (u'(x),v'(x)) are the 
solution of the following rational interpolation problem: 

with 

f U n = 1, V n = 0, 
1 < = °> V '„ = !• 

Definition 4.5. The r-degree of a vector polynomial w(x) = (wi(x) ui2(x)) T is defined as 

t — degw(x) :— max{degu>i(a;), degu^Oc) — t} 

Definition 4.6. A polynomial vector in K[x] 2 is called r-reduced if the r-highest degree coefficients 
are linearly independent. 

By construction, the columns of B(x) form a n-reduced basis of the module of polynomial vectors 
r(x) G K.[x] 2 that satisfy the interpolation conditions 

/ir(wj)=0, j = 0,...,2n-l 

with fj = (T(ujj),uj™) G K 2 . The columns of B(x) are also called a n-reduced basis for the interpolation 
data (u>j,fj), j = 0,...,2n- 1. 

Theorem 4.7. Let r = n and J be a positive integer. Let Ai, . . . , Aj G K and <f>\, . . . , cf>j G K 2 \ 
{(0,0)}. Let 1 < j < J and tj G Z. Suppose that Bj(x) G K[x] 2x2 is a Tj-reduced basis matrix 
with basis vectors having tj — degree Si and 62, respectively, corresponding to the interpolation data 
{(Ai,0i);i = 

Let Tj^j := 61 — 62- Let Bj^j(x) be a Tj-yj-reduced basis matrix corresponding to the interpolation 
data {{\i,(f>i Bj(Xj));i = j + 1, . . . , J}. 

Then Bj(x) := Bj(x)Bj^.j{x) is a tj -reduced basis matrix corresponding to the interpolation data 
{(Xi,(t>i);i = 1, . . . , J}. 
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Proof. For the proof, see [20] , 



□ 



When we apply this theorem with Xj — ujj £ i^n as interpolation points, we obtain a superfast 
algorithm in O(n\og 2 n) to compute B{x). See [5D] for more details. 



5. Conclusion 

In this paper, we re-investigate the solution of a Toeplitz system T u = g from a new point of view, 
by correlating the solution of such a problem with generators of the syzygy module C(T(x) , x n , x 2n — 1) 
associated to the Toeplitz matrix T. We show that C(T(x), x n , x 2n — 1) is free of rank 2 and that 
it has a n-basis. We show that finding a n-basis of C(T(x), x n , x 2n — 1) is equivalent to computing 
the Gohberg-Semencul decomposition of T _1 , and we reduce the solution of T u — g to an Euclidean 
division. We give two superfast algorithms computing a n-basis of C(T (x) , x n , x 2n — 1) and a superfast 
algorithm to obtain the solution from this n-basis. 

A perspective of this work is to generalize the approach to two-level Toeplitz systems or to Toeplitz- 
block- Toeplitz matrices and to correlate the basis computation of a multivariate syzygy module to 
"Gohberg-Semencul" decompositions for Toeplitz-block- Toeplitz matrices. 
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